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Abstract 

In this paper the existence and unicity of a stable periodic orbit is proven, for a class 
of piecewise affine differential equations in dimension 3 or more, provided their interaction 
structure is a negative feedback loop. It is also shown that the same systems converge toward 
a unique stable equilibrium point in dimension 2. This extends a theorem of Snoussi, which 
showed the existence of these orbits only. The considered class of equations is usually studied as 
a model of gene regulatory networks. It is not assumed that all decay rates are identical, which 
is biologically irrelevant, but has been done in the vast majority of previous studies. Our work 
relies on classical results about fixed points of monotone, concave operators acting on positive 
variables. Moreover, the used techniques are very likely to apply in more general contexts, 
opening directions for future work. 

1 Preliminaries 
1.1 Previous works 

Among recent biological concepts, gene regulatory networks is a particulary intriguing one. 
Coarsely said, it refers to a set of genes coding for proteins being able to activate or inhibit the 
expression of other genes within the same set. Since such systems are likely to involve large 
numbers of genes, and since moreover gene expression is known to be a non- linear mechanism, 
mathematical models are, though a required tool, facing a major challenge. This fact has 
impelled a large amount of work, giving birth to numerous models, see [2] for a review. Among 
different types of models, a special class of piecewise linear differential equations have been 
proposed in the early 1970's [S]. This class has the advantage of being mathematically much more 
tractable than more classical, smooth non-linear equations, nonetheless offering a comparable 
variety of behaviours. 

Since its introduction, the class of picccwise-lincar models has thus been studied in various 
fashions. Periodic trajectories, or limit cycles, of these systems have been especially investigated. 
The first results concerned a particular case, characterized by the equality of all decay rates of 
the system. Under this assumption, trajectories are straight lines in regions where the system 
is affine, and the whole analysis becomes much simpler. Notably, an explicit Poincare map is 
possibly derived, which is a fractional linear map. It is then possible to reduce the problem 
of existence and stability of periodic orbits to an eigenanalysis. This method was introduced 
in [SI [7] , and subsequently improved by different authors [TH [HI [S] ■ 

To our knowledge, ony one study did not require this uniform decay rates assumption, but 
concerned on the other hand systems whose interaction structure consists in a single, negative 
loop [14] . Biological examples, such as the repressilator [3] make this type of interaction structure 
relevant. For such systems, the existence, but not unicity, of a stable limit cycle in dimension 
^ 3, and of a stable focus in dimension 2 is proven in |14j . This proof relies on Brouwcr's fixed 
point theorem, and is not entirely self contained. 
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Our main result is a complete proof of the existence and unicity of the mentioned attractors 
in feedback loop systems. Moreover, the mathematical tools involved to prove this result seem 
to have been ignored in the context of piecewise linear gene network models, and offer a very 
promising framework. 

More precisely, the result on which we mainly rely is a fixed point theorem for monotone, concave 
operators, where monotonicity and concavity are defined with respect to a partial order. There 
exists actually a wcU-developped corpus about operators satisfying such properties, especially 
for systems with positive variables. These results have the advantage of being generic, and 
the hypotheses on which they rely are on the other hand naturally satisfied in the context of 
piecewise affine models of gene networks. In order to give a more precise foretaste of the rest of 
this paper let us state now the main result. To make it understandable, let us just say that C is 
a periodic sequence of regions, crossed by the trajectories of any piecewise affine model of gene 
networks consisting in a negative feedback loop. Also, is a particular subset of one region, 
or 'zone', which is a well defined Poincare section, and T is a first return map on Z^. 

Theorem 2. Let T : Z^ ^ Z^ be the Poincare map of 6. 

• If n = 2, then V.t G Z^, T™x when m ^ oo. 

• If n > 2, then there exists a unique nonzero fixed point q — Tg. 
Moreover, q^Z} and for every x £ Z^ \ {0}, T™x — > g as m ^ oo. 

A rapid overview of the tools and references we rely on is given in section [TT^ Then, the 
class of piecewise affine models we focus on is presented in section [51 Some general observations 
about monotonicity and concavity for this class of systems are made in section[31 These remarks 
lead to focus on negative feedback loop systems in section [H where the main results are stated 
and proven. 



1.2 Useful mathematical results 

The main result of this paper strongly rely on one theorem, that we recall here. It is a fixed point 
theorem of monotone, concave operators. Both monotonicity and concavity shall be understood 
with respect to some partial order, see below for definitions and notations. Operators with both 
these properties have lead to various results on existence and unicity of fixed points. Among 
early studies, Krasnosel'skii has provided important results [TT]. Here, we use a more recent 
presentation of such results, due to Smith [13], whose formulation is well adapted to our purpose 
and context. Namely, we use theorem 2.2 of [T5] . 

Before stating the latter, let us introduce some notations (which are similar to those in [T5]). 
We denote x < y and x ^ y if these inequalities hold for each coordinate (resp. entry) of 
vectors (resp. matrices) x and y. Then, we denote x < y ii x ^ y and x ^ y. For x ^ y, 
[x,y] = {z \ X ^ z ^ y}, and {x,y) = {z\x < z < y}. For any set A, A denotes the interior of 
A, and ci{A) its closure. Now, we can state the following, 

Theorem 1. Let p e R" , and T : [0,p] — > [0,p] continuous, in (0,p). 
Suppose DT{0) = lim DT{x) exists. Assume: 

x>0 

(M) DT{x) > if X > 0, X <p. 
(C) DT{y) < DT{x) if < x < y < p. 
Assume also Tp < p. 

Suppose that TO = 0, and define A = p{DT{Q)), the spectral radius of DTifi). Then, 
A ^ I V;r £ [0,p], r".T when ri — > oo. 

A > I There exists a unique nonzero fixed point q = Tq. Moreover, q G (0,p) 

and for every x £ [0,p] \ {0}, T"x — > g as n ^ oo. 
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The case TO ^ leads to either a unique fixed point, or to diverging orbits. However, this 
case will not appear in our context, and thus it is not detailed in the theorem above. 
One may remark now that in the case when T is twice differentiable, the concavity condition 
(C) admits a simple sufBcient condition. The map Ti : [0,p\ — > R+ denotes the ith coordinate 
function of T : [0,p] ^ [0,p]. 

Proposition 1. Suppose that for all i, j, fc £ {1 • ■ • n}, and for all < x < p, 

dxkdxj ^' 

and for all i^j there exists a k such that the inequality is strict. 
Then T satisfies condition (C) of theorem{I\ 

Actually, it is clear that under this condition each term of DT is a decreasing function 
of each coordinate Xk- It is moreover strictly decreasing in at least one of these coordinates, 
and (C) thus follows. Observe by the way that the notion of concavity (w.r.t a partial order) 
we deal with here is weakened by the fact that it concerns only ordered pairs {x, y) of variables. 
Remarkably, theorem [T] relies on Perron- Frobenius theorem, which will be used at some point 
in the proof of our main result. It usually says that any matrix with positive entries admits 
a positive, simple eigenvalue, associated to a positive eigenvector. Here, we will only need the 
following corollary: if A is a real n x n matrix such that A > 0, then A admits a positive real 
eigenvalue. 



2 Piecewise afRne models 
2.1 Formulation 

The general form of the piecewise affine models we consider here may be written as: 

^ = n{x)-Tx (1) 

The variables {xi . . represent concentrations in proteins or mRNA produced from n inter- 
acting genes. Since gene transcriptional regulation is widely supposed to follow a steep sigmoid 
law, it has been suggested that idealized, discontinuous switches may be used instead to model 
these complex systems [5]. Accordingly, s+(- ,6') : R ^ {0, 1} denote the increasing step func- 
tion, or Hcavisidc function: 

s+{x,e) = if x < 6*, 

s+(.T,6i) = 1 ifx^e, 

which represent an effect of activation. Also, s~{x,6) = 1 — s'^{x,9), is its decreasing version, 
and represents inhibition. 

Then, k : R" R" is a piecewise constant production term that can be expressed in terms of 
step functions s^{xi,6i). T E R"^" is a diagonal matrix whose diagonal entries Tu = 7^, are 
degradation rates of variables in the system. 

Since each variable Xi is a concentration (of mRNA or of protein), it ranges in some interval of 
nonnegative values denoted [0, max,]. When this concentration reaches a threshold value, some 
other gene in the network, say gene number j, is suddenly produced with a different production 
rate : the value of kj changes. For each i € {1 ■ ■ - n} there is thus a finite set of threshold values, 
which serve as parameter of step functions : 

e,^{el...,er'}, (2) 

where the thresholds are ordered: < 9} < ■ ■ ■ < B^^'^ < maxi. The extreme values and max^ 
are not thresholds, since they bound the values of Xi, and thus may not be crossed. However, a 
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convential notation will be : 0° = 0, and 91' = maxi. 

Now, at a time t such that Xi{t) € Qi, there is some j G {l---n} such that Kj{x{t'^)) ^ 
Kj(x{t~)). It follows that each axis of the state space will be usefully partitionned into open 
segments between thresholds. Since the extreme values will not be crossed by the flow (see 
later), the first and last segments include one of their endpoints : 

^^ e {[el el), {er\ ef]\ u {(^^ ei^')\j e {i • • • - 2}} u (3) 

Each product T) ~ JlILi defines a rectangular domain, whose dimension is the number of 2)^ 
that are not singletons. When dim D = n, one usually says that it is a regulatory domain, or 
regular domain, and those domain with lower dimension are called switching domains [3j. We 
use the notation ^ to represent the set of all domains of the form above. Then, &r will denote 
the set of all regulatory domains, and &s the set of all switching domains. 
A convenient way to represent regular domains makes use of a discrete map, whose range we 
write as: 

n 

A = \{{Q---q,-\}. (4) 

i=l 

The map in question is then defined as d : — *■ -A, and sends each regular domain to the 
superscript of its "lower-left" corner: d ij\i{el'~^ , 0"')) = (ai — 1 . . . a„ — 1). In the following, 
we will often identify regular domains and their image in A, talking for example about some 
'domain a'. 

The dynamics on regular domains, called here regular dynamics can be defined quite simply, 
due to the simple expression of the flow in each D € On sets of Sig on the other hand, the 
flow is in general not uniquely defined. It is anyway possible to define solutions in a rigorous 
way, yielding what will be mentioned as the singular dynamics, at the price of considering set- 
valued solutions. The latter's definition rests on Filippov's theory of differential equations with 
a discontinuous right-hand side, and its formulation in the particular case of equations U]) can 
be found in [T1|TU], among others. Since we will only encounter systems for which the Filippov 
solutions lead to single valued trajectories, we shall only describe the regular dynamics with 
some detail. 



2.2 Regular dynamics 

Regulatory domains are of particular importance. They form the main part of state space, 
and the dynamics on them can be expressed quite simply. Actually, on such a domain D, the 
production rate k is constant, and thus equation ([T]) is affine. Its solution is explicitly known, 
for each coordinate i : 

ip^{x,t) = a;j(t) = — - e"''-* [xi- — ], (5) 
and is valid for all t E R-|_ such that x{t) E T). It follows immediately that 

<A(2)) = (</>i...<^„)=f^...^) 

V 71 In J 

is an attractive equilibrium point for the flow Hence, if it lies inside D, it is an actual 
equilibrium of system ([T]). Otherwise, the flow will reach the boundary dT) in finite time, unless 
0(D) lies exactly on the boundary of D. However this situtation is clearly not generic, and we 
assume in the rest of this paper it never occurs. At the time when the flow reaches dT) thus, the 
value of K will change, and that of (f) accordingly. The point (/'(D) is often called focal point of 
the domain D. Then, the continuous flow can be reduced to a discrete-time dynamical system, 
with a state space supported by the boundaries of boxes in 

Hence, the state space of this discrete-time system is part of i^,,, which may seem problematic 
at first sight. In fact, it is actually not always possible to define discrete-time trajectories on 
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domains of having dimension n — 2 or less. On any n — 1 domain D, on the other hand, such a 
definition may always be properly provided, in terms of the flow lines in the two regular domains 
separated by D. If these flow lines both point towards, or away from D (which is expressed 
formally using the flow coordinate which is normal to D), the latter is called respectively black 
wall or white wall. Otherwise, i.e. when flow lines both cross D in the same direction, one 
usually refers to 2) as a transparent wall. In the two flrst cases, the Filippov theory is required, 
providing set-valued trajectories on the walls, the so-called sliding modes jlOj. On transparent 
walls, trajectories are simply deflned by continuity, from the flow lines on both sides. We shall 
encounter only this third kind of wall in the rest of this paper. 

Note that if D G is represented by a G ^1, i.e. d{D) = a, one shall most often denote by 
4>{a), or 0", the focal point associated to this domain. 

Once the flow ^ is given in a box T)a, it is easy to compute the time and position at which 
it intersects the boundary of Dq, if ever. The possibility for each facet to be encountered by the 
flow depends uniquely on the position of the focal point : {x\xi ~ ^'f'^^} (resp. {x\xi ~ 6*°'}) 
can be crossed if and only if < (resp. 0i > 0'^'). According to this observation, we 

denote /+,(a) = {i € {1 ■ ■ -njl 0, > ^'j, and - {* G {1 • • -ri}! 0. < ef-'}. 

Then, Iout{o) ~ I'^^^{a) U I^^^{a) is the set of escaping directions of Dq. 

Also, it will be practical to point out the lower and upper thresholds bounding a box, thanks 
to the pairs of functions Of : yi — > O^, 9~ {a) = 61°'"^ and 9^ {a) = 61\ These functions will 
clearly only be useful notations, and bring no further information about the system. 
When it is unambiguous, we will omit the dependence on a in the sequel. 

Now, in each direction i G lout the time at which x(t) encounters the corresponding hyperplane, 
for X G Da, can easily be shown to be: 

r.(.) = ^lnfmin(^Vffi,^i^^ 

li \ I (pi- Xi (pi- X, 

Taking the minimum t(x) = minjg/^^j Ti{x), and reinjecting it in equation (0), we get the exiting 
point of Da when starting at x. Since this process is intended to be repeated along trajectories, 
X will generally lie on the boundary of the current box, except for the initial condition, which 
may however be chosen on a facet without loss of generality. We then get a transition map 
7'' : dDa ^ dDa : 

J^X = ip{x,T{x)) , . 

= (j) + a{x){x^(t)). 
where a{x) = exp(— r(x)r). The latter depends on a, as seen from ([6]). 

Now, the initial system ([1]) may been reduced to a discrete time dynamical system, as mentioned 
earlier. It consists in iterates of a global map T on a domain Domo". As already mentioned, 
Dom ^ must be contained in the set of n — 1 domains of . Since T has to be iterated on this 
domain, all points which reach some n — 2 dimensional (or less) domain after applying T a finite 
number of times, must be excluded from Domg-. Furthermore, whereas T is readily defined on 
transparent walls, black and white walls are much less obvious to deal with. We thus suppose 
from now on that Dom^j entirely consists of transparent walls. A more detailed discussion about 
the topology of this domain can be found in [5]. A simple criterion to ensure that all walls are 
transparent is the absence of autoregulation, in the sense that no production term Ki depends 
on Xi. 

At each point of such a domain, it is possible to consistently apply a unique local map T". 
Actually, each wall is either on the boundary of the whole domain IJ^ Da, and thus of a single 
regular domain Da, or can be written as an intersection dDa H dDi, of two regular domain 
boundaries. It is clear that on transparent walls there is always exactly one of the two domains, 
say Da, such that T" is not the identity. Applying T° is then consistent with the orientation 
of the flow lines. Yet, we only care here with forward trajectories, since is not properly 
defined on the full domain Domg-. 

Suppose now that there is a wall W C Domg-, and a sequence . . . of regular domains such 
that (T°' o T"'~' • ■ • o 7°-^){W) nW ^ 0. Now, it is tempting to study the dynamics of this 
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iterated map on W . In particular, any fixed point of this map determines a periodic trajectory 
of the continuous time system ([T]). 

3 Regions of fixed monotonicity and concavity 

In order to apply theorem [T] to an iterate of the transition map, as defined in previous section, 
it is necessary to check whether it is monotone and concave, in the sense of properties (Af) and 
(C). In this section it is shown that these conditions cannot be satisfied by the transition map 

on a whole wall, and that some refinement should be done. To this aim, one should compute the 

I i—i 1 

jacobian of a composite map, of the form: £'(T'' oT° •••oT° ). This requires the computation 
of an arbitrary DT" . 

Let W be an arbitrary wall. One has seen in previous section that T = T" , for some uniquely 
defined a. Let moreover x & W such that there exists a unique s € Iout{a.) such that t{x) = 
Ts{x), and denote W C {x\xs — dg} the wall such that 7x G W . Then there exists an open 
neighbourhood of x in which the following holds: 

7,x==(bf + ix,~(j)f)a,{x), (8) 

where is the ith coordinate function of T. To abbreviate computations, we have denoted 
the ith diagonal entry of the diagonal matrix a{x) defined in previous section: 



a.i{x) = 




Then it is rather straightforward to compute the partial derivatives at a point x: 

{{asix))"^ if i=j 

-/l (t)f -Xi 2^ . (IQ) 

otherwise. 

It appears from pO|) that the partial derivative are of constant sign if, and only if, sign {(pf — Xi) 
is fixed. Actually, all degradation rates 7.; are positive, and as{x) = exp(— 7sr(x)) S (0, 1) since 
t{x) is nonnegative. 

Now make two observations on the sign pattern ol (jf" — x in general. The first one concerns 
its invariance with respect to time evolution: 

Proposition 2. From the expression 0), one can observe that the quantities sign (0° — Xi) are 
preserved by the flow, in the sense that sign [(pl — Xi) = sign (0° — (pi{x, t)) for any t ^ 0. 
In particular, since the wall W chosen at the beginning of this section is included in a hyperplane 
of the form {x \ Xj = 6j}, the above observation yields sign (0" — Tjx) = sign (0° — 9j). 
Similarly, for any x which escapes via the wall W C {x | Xs = 0^}, one has sign ((/)° ~ Xg) = 
sign{rs-Os). 

A second fact concerns the influence of the different coordinates on sign ((/)° — x). 

Proposition 3. For any escaping direction i e Iout{o-), one has either Xi < {a) < 0°, or 
(jjf < 0~ (a) < Xi, and in both cases sign {(jjf — Xi) is constant in the whole closure c£(T>a) of the 
domain. IfOi e {d^{a), 9^{a)} is the escaping threshlod, this sign equals sign(0° — Oi). 
For the non escaping directions on the other hand, sign ((/)° — Xi) is not constant on W , since 

e-[a)<<t>'} <et{a). 
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Thus, condition (M) cannot be valid on the whole wall W, unless all directions are escaping. 



Since a differential of the form Z?(T° o T° • • • o ) is a product of differentials of local 
maps D{7°' ). and since the latter may only be positive in a region where the vector 0" — x has 
a fixed sign pattern, it seems relevant to restrict the study to such regions. 
Now, condition (C) concerns the variations of the terms in pO[) . As shown by proposition [1] 
the second-order derivatives are informative in this case. They can be expressed as: 



7j {as{x))-<' 



Is 

2i 

Is 





if 



J, m 



or J = s, TO 



{as{x))''^ if m=j = s 
otherwise. 



(11) 



Here again, these derivatives are of fixed sign only if 0" — x has a fixed sign pattern. Hence, 
both first order and second order derivatives are of constant sign in the same regions. One is 
thus naturally led to consider rectangular regions partitioning the wall W^, defined by the sign 
oi (j)"" — X, or in other words by the position (in terms of the partial order) of points with respect 
to the focal point. 



It will be shown in the following sections that in the special case of negative feedback loop 
systems, only one such region contains the long-term dynamics, on each wall of a periodic 
sequence of domains in state space. Before proving this, let us discuss a slightly more general 
case than feedback loop systems. Our hope here is that this discussion give some intuition of 
the basic facts that will serve afterwards. It also suggests that tools from positive operators 
theory might be efficient in a broader context than the sole feedback loop systems. 
Namely, suppose given a periodic sequence of domains , . . . a^,a^~^^ = , such that there 
exists a unique escaping direction for each box a*. This case is actually the simplest one 
when studying periodic orbits, since if there are several escaping directions for one box in 
the sequence, then some trajectories escape the cycle at this box. The set of points remaining in 
the periodic sequence . . .a^ is not easily defined in general, though under the uniform decay 
rates assumption it may be expressed by a list of affine inequalities [U E] • 

Then, each wall reached by successive iterates of the transition map has to be partitioned into 
2n-2 j-QgjQjjg jj^ general. Actually, on such a wall points have n — 1 coordinates which are not 
constant: in one direction it equals a threshold value. Among these n — 1, one exactly is the 
escaping direction Si_i of the previous box, a'^^. Thus according to proposition [21 the following 
is fixed 

sign - J = sign - J, 

for some to, and any point x on a trajectory crossing the considered sequence of domains. 
Thus if one assumes Si ^ Si-i, then n — 2 directions remain for which two signs are possible. 
At each wall, the images of these 2"~^ region may intersect several regions among the 2"~^ 
partitioning the next wall. Hence, one may refine the partition of the latter, and iterating this 
process, obtain a partition of a wall into regions such that all partial derivatives of iterates, of 

the form -g^, for fc G {1 ■ ■ ■ 1} are of fixed sign, and the second order derivatives as well. This 
process would lead to consider a high number of regions in general, each of which may satisfy 
the hypotheses of [T] This seems far from trivial in the most general case, but would certainly 
be fruitful in some more particular cases. 

Hence, let us look for such special cases of equations ([T]). For two successive boxes a*~^ and a*, 
one has, as seen from ((HJ: 

0f - 7,x = 0f - ^f" + i^f" ~ X,) a, (x), (12) 

and if the focal points are 'aligned', in the sense that 4'j=4''j : the sign of the quantity 
4>j — Xj remains unchanged after applying T. Actually one as seen already that as{x) g (0, 1), 
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and applying T implies that the focal point must be replaced by that of the following box. 
This geometric condition should thus avoid us to refine the partitions of successive walls, as 
mentioned in the above discussion, since it imposes that these regions do not overlap from one 
wall to the following. 

It may be stated more formally. Recall that a* refers to the ith crossed box, and that Si is its 
unique escaping direction. Assume a strong version of the previous alignment of focal points: 

Vze{l---£},Vje{l---n}\{s,+i}, cj/f = c^^' ■ (13) 

This alignment assumption will allow us to use theorem [1] in the region where the long-term 
dynamics take place. Furthermore, it will be shown that the geometric property (|13p can also 
be obtained as a consequence of further assumptions on the interaction structure of the system: 
if the latter consists in a single negative feedback loop, it follows that there exists a periodic 
sequence of boxes satisfying (fT^ , as we will see in the next section. 



4 Negative feedback loop systems 

4.1 Preliminary definitions 

In this section, we focus on systems of the form: 

(14) 

K° + k} S+(.T,;_i, 6'|_ J - JiXi, 1 = 2. ..n. 

In such systems, each variable Xi activates the production of the next variable, Xi+i, except Xn 
which inhibits xi. Hence, their interaction structure is that of a negative feedback loop. Up to 
a translation, they are of the form dealt with by Snoussi in [T3]. In the case n = 3 one finds the 
well-known repressilator [3]. Remark now that any production term Ki is independent on the 
variable Xi, and thus all walls must be transparent. 

Without introducing unnecessary formalism, let us only say that we call interaction graph, IG is 
an oriented, labeled graph, whose vertices represent variables of the system, and edges represent 
interactions between them. For example, Xi — >~ Xj if an increase of Xi may lead to a decrease 
of Kj at least at some point x in state space. Then, a negative loop is a cycle in IG with an odd 
number of negative edges. It is easily shown that any system whose interaction graph consists 
in a single negative loop (involving all variables) is equivalent to up to a symmetry in state 
space. 

The particular form of (fT4| implies that there is only one threshold 6} in each direction i, and 
0°, 9f represent the bounds of the range of Xi, in accordance with 
In the following we always assume that 

V^ e {1 • • • n} K° + kI > 7, 0} and k"^ < 7, (15) 

This implies that higher values of focal points are above thresholds, and lower values are below 
thresholds. Otherwise, all trajectories would converge to a unique fixed point, see [14j lemma 1. 

Before entering into more detail, let us operate a few simplifying operations. First, we chose 
to set Wj = Xj — 6*1 . Then, 

Wj = n° + K]h=^(wj_i, 0) - 7j6'] - 7jWj, 

and thus all trcshold values are zero. Accordingly, all regular domains arc now identified with 
(bounded rectangular regions of) the 2" orthants of R". They arc mapped as explained in 
section \TT[ onto the discrete set A = {0, 1}". The corresponding discrete mapping d : — > -A 



' dxi 

IT 

dxi 

- 'dt 
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sends negative coordinates to 0, and positive ones to 1. 
The bounds of the whole domain arc now denoted 

Oi^e^-ol^-el and e+^el-el 

Observe that variables may now be negative, and do represent concentrations only up to a 
translation. 

The new dynamical system behaves identically to p^ . once thresholds have been set to zero, 
and the following renamings have been done: 

^ K° — 7j;0i , and ^ </>,; — 9}. 

Then, for each i £ the focal point coordinate 0i may only take two values. We 

abbreviate them as follows: 

4^- = ^ and </)+ = ^!l±^. (16) 

li li 

One shall also note = 7^ 0^ . 

Remark that 0~ is negative and 0^ positive (and similarly for k^). From this observation, one 
deduces that the focal points are entirely determined by the sign of their coordinates. In other 
words, the coordinates of any focal point are known exactly from the regular domain it belongs 
to. More precisely, if a focal points belongs to a domain a <^ A, the sign of 0i is lai — 1 e {±1}, 
for each i £ {1 ■ • • n}. 

Now it is time to describe the dynamics of . In brief, we now focus on systems of the 
form below, into which any system (|14p can be transformed: 

+ (k+ - Kj;)s"(a;„,0) - 71X1 

(17) 

K- + (k+ - K- ) s+(xi_i, 0) - 71X1, i = 2...n. 

A first, coarse-grained description can be achieved in terms of discrete transitions between 
qualitative states in A. This can be conveniently explained with the aid of the state transition 
graph, or simply transition graph, whose vertices are elements of A, and edges (a, b) exist if and 
only if there exists some continuous trajectory crossing successively the corresponding regular 
domains in state space, Dq and Db. This graph will be denoted TG. 

It is easily deduced from the properties presented in section [^T^ that TG is entirely determined 
by the position of focal points. More precisely, (a, 6) S £ if and only if 6 = a = <p{a), or 

b E {a + Bi \ i such that = and 0i(a) > 0} U {a — | i such that = 1 and (t>i(a) < 0}. 

Observe that a and 6 may differ at one coordinate at most, since TG represents the regular 
dynamics only. 

From (|17p and p3)) . one can deduce the presence of the following cycle of length 2n in TG: 
1---11 — > 01---1 — > ■•■ — > 0---011 — > 0---01 

I 1 

1---10 < — ••■ < — 110 •••0 < — 10- --O < — 00- --o 

It can be shown that there is no fixed state in TG, and that C is the only cycle without escaping 
edge. Then, we denote 





a2 






a" 








a2" 


1---11 


01 •••1 


001 •••1 




0---01 


0---0 




1 ■ • • 100 


1---10 



(18) 



Furthermore, let us write (jf = (/'(a*), the focal point of domain a\ and the wall between 
boxes a' and a'+i, i.e. it verifies c£{W') = cl{Da^) n c£(2)a.+i). 



' dxi 

~dt' 
dxi 

. ~dt 
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Figure 1: The transition graph of a feedback loop system of the form 07)) with three 
variables. The cycle C is represented with plain lines, and the walls W" have been 
schematically depicted as well. 



In the following, the superscript in a* must be understood periodically, so that for instance 
0^"+^ will mean a^. It appears that each pair of successive states a', a*"*"^ in this cycle differ 
only at one coordinate, which we denote Si. If moreover we denote sf where ± is the sign of 
al'^^ — al., which is also the sign of 0*., one obtains 



i 


1 ■ 


n 


71+1 • 


2n 




1- ■ 


nT 


1+ . 


■ n+ 



(19) 



Hence Si is the only escaping direction of box a', a fact that can also be written as loutW) ~ {si}- 
As seen in section r2.2[ this implies that all trajectories in a domain a' leave it in finite time, via 
the wall VF* C {x | = 0}, and enter domain a*+^. It follows that no trajectory can escape 
the cycle, i.e. C is an invariant region. This cycle is represented on figure [1] for n ~ 3. 



A rapid inspection shows that the walls can be described in a quite explicit manner: 

n 0) X {0} X n (0, 9+] for * e {1 . • • n}, 

W'^{ ^ (20) 

n ] X {0} X n [^7' ^) for i e {n + 1 • • ■ 2n}. 

j<Si j>Si 



From the remark under equation ()16p it is not difficult to deduce that each pair of focal 
points 0', (f)^'^'^ differ only at coordinate s^+i. Actually, they verify G D^i+i, and one has 
seen that a'+^ and a'+^ only differ at coordinate s^+i. Hence, the cycle 6 satisfies the alignment 
condition (|13p seen in previous section. We can state it with our current notations: 

V*e{l---2n},Vje{l---n}\{s,+i}, 0} - 0}+^ (21) 

Remark also that by construction one has for the remaining coordinate: 

V*e{l..-2n}, 0:,^, (22) 

Since our aim is to study limit cycles arising from p4p . with the help of theorem [I] it seems 
natural to consider the following map: 

T°''J"'"T"'"~' • ■ - T"' : — >W^ 

Although each T" is defined on the whole boundary of T)a, see (O, one considers here the 
restrictions to walls W^~^ only, which we write explicitly: 

ja" : iy*-i — > W 

(23) 



j=l...n 
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where, as in ((H]), 

Remark here that the expression of aj depends on the number i of the region under consid- 
eration. However, this region wiU always be clear from the context, and we thus leave this 
dependence hidden, avoiding notations such as a* (a;). 

4.2 Restriction of domains 

4.2.1 Projection and extension to the closure of walls 

To fit the hypotheses of theorem[Tl one must consider a closed domain of the form [0,p], p E R" , 
with nonempty interior, and since the walls are n — 1 dimensional, their interior is empty. 
Therefore, one may first project to R""-'^, and drop out the first coordinate, which is zero. 
Similarly, there is no loss in suppressing the Sith coordinate of points in W'\ since it equals zero 
by construction. 

Second, each has an extension to the closure c^(W^'~^), defined by continuity of the expres- 
sion To compose these extensions, one has to check that each of them maps c^(VF*~-'^) to 
c^(VF*). Since the walls themselves are correctly mapped, and are products of one singleton {0} 
and of segments of the form (0,0^] or [9~ ,0), it suffices to check that the points of c£(W^~^) 
having at least one coordinate equal to are mapped to (£(W^). 

So, let X G ci(W^~^) such that Xj = 0. Since c£(W^~^) C {x \ Xsi_i = 0} it is relevant to 
consider only j ^ s^-i. Then, T" x g ciiW"^) if, and only if, each coordinate x belongs to 
the projection in direction k of diW"^), which is either [0, 6*^] or [0^, 0]. Let us call this segment 
"Jk to simplify the exposition. The explicit formula given in ()23p shows that the coordinate map 
Tj depends only on variables xj and Xs^- Hence, if j ^ Si, one only has to check whether 
7j X E 3j , while if j ~ Si all coordinates must be checked. These two cases give 

• li j Si, T^'x = (1 — C(j{x)) (f)j. From the inequalities between Xj, 0* and 9j = 0, it is 
easily checked that aj{x) £ (0,1]. Moreover, for j ^ {si} = Iout{a'^), 4>] is in 3j \ {0}. 
Then, 7fx e follows easily. 

• If j = Sj, remark first that x G d{W''-^) n c£{W') . Then 7f^x = G 3j. Moreover, for 
any fc G {1 ■ ■ • n}, as^. (x) = 1, so that x — x^, that is to say a; is a fixed point. 

It follows that in any case, T'^ x G cf (W^*), as awaited, and we consider now that the maps T" 
are defined on the closure of walls. 

Finally, proposition [5] implies that the image T" (cf(W^'~^)) is always contained in the subset 
of ci(W^) whose points satisfy 

sign -Xs^_,)=sign{<f>l^_J. (24) 

This restriction is illustrated on figure [21 Since our interest is in the iterates of those maps, we 
should restrict them to such subsets. 
In summary, three simplifications can be done: 

• The Sith coordinate of x G VF* can be suppressed without loss of information 

• The transition maps can be defined on the closure of walls. 

• Only the points verifying ([24]) can be reached by . 

The notation will be used to denote the domains obtained after these three simplifications. 
More explicitly is obtained from (|20p by suppressing the {0} term, using closed intervals. 
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Figure 2: The same walls as in figure [TJ this time in phase space. Dotted lines 
represent thresholds. The dashed lines on a wall represent lines where Xj = (pj. 
The alignment condition if^) is fulfilled on the left, and not on the right hand side. 
The right-hand side figure is thus only presented as a counterexample to lUp. The 
blue region on each such wall represents the image T° (W'^~^). They are shown 
as polyhedras for ease of observation, but they have in general a piecewise smooth 
boundary, and not a piecewise linear one (unless all decay rates are equal). The 
shaded regions are not reachable, as an illustration of eq. |23)- 



and replacing moreover [dg.__^,0] by [0s;_j,O], or [0, 6'+_J by [0, 0+_J. For example: 

n n— 1 

= {0} X J](0, e+], so that W^' - n ^ ['^' '^"]- 

1=2 i=2 

Each is a rectangular region in R"~^ with nonempty interior, still called wall in the 
sequel, and the preceding discussion shows that each map T° induces a well defined map 

Remark that although points in W^~^ do not have a Si_ith coordinate, is defined using 

expression ()23p . with Xsi_-^ — 0. Let us write it explicitly for later purpose: 

?1,^= (1 (25) 

Now, we define a map 

T = T'''t°'"T"'""' • • • : — >W^ (26) 
Our aim will be to check the hypotheses of theorem [T] for this map T. 



4.2.2 Partition of walls 

Let us now consider how may be partitioned according to the sign pattern of (/>* — x. 
Following (pi)) . the coordinate Si_i of this sign is fixed by construction. As for the coordinate 
Si, it has been suppressed from by construction. On the other hand, proposition [3] implies 
that this sign is not constant for any other coordinate j. Then, denoting 

J; = {1 • • - n} \ {s.t_i, sj, 

it follows that the above mentioned regions may be described using sign vectors in the finite set 
= {il}"^'. Vectors in will then have their coordinates indexed by Ji, which will prove 
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more convenient than indexing them by {1 • • • n — 2}. For any sign pattern cr G E^, let us define 
a corresponding 'zone' 

r{a) = {x e I Vj G J^ , sign (0j - x,) - ^r,}. 

Now, from the translation of all thresholds to zero, the origin of R"~^ belongs to every wall 
W^. Moreover, no focal point have a coordinate equal to zero. Hence, on each wall W^, there is 
a unique zone containing the origin on its boundary, and the corresponding sign vector in is 

= (sign(0}))^^^^ (27) 

We will see now that these particular zones, having the origin on their boundary, in fact attract 
all trajectories of the system, and then remain invariant. As they will be the only zones we 
consider eventually, we simply denote: 

= Z\a'). 

Remark that the definition of these zones and equation (|20p allow for a more explicit formulation: 

Z' = \ ^ (28) 

n '^j ] >^ n ["^^J^' 0] fori G {n + l---2n}. 

I j<Si j>Si 



4.2.3 Attractive and invariant regions 

Proposition 4. The zone Z^ is attractive for the dynamics induced by T: 

T (w^'j C Z\ 

Proof. From ^ and dH]) it follows that sign {(pi) = (-1, 1 • • • 1) G {±1}". Then, si = 1 and 
So — 32,1 = n implies Ji = {2 • • • n — 1}. Then, tr^ = (1 • • • 1) G {±l}-^\ and 

Z^ = {x G I Vj G Ji , Xj < 0]} 

71. 

= nt'^' since 0] = ^+ for j G Ji. 

Let x^G W^. In the following we denote = T°'a;'~^ for i ^ 2, so that Tx^ = Clearly, 
x' G W for all i. 

Now, for each i G {1 • • • 2n}, the alignment condition ([2T|) and the expression ((23|) of imply 
the following conservation of signs: 

G {1 • ■ • 2n}, V? ^ s,+i, sign (0j - xj) = sign (0}+^ - x}+i). 

In particular for the n last domains of the cycle 6, remarking from table p9p that s„+i = i, the 
above has the following consequence: 

Vi G {1 • • • n}, V? <i + l, sign (?i}'+' - x]+'') = sign (0;'+'+^ - x]+'+'^). 

This equality can be propagated, in the sense that 

sign (0^1 - X^l) = sign (0^2 _ ^«+2^ ^ ^ ^ ^ ^ gig^ (^2n „ ^2„) ^ gig^ (^1 „ ^2„+l)^ 

sign (</>^+2 _ ^ gjg^ _ ^„+3) ^ ^ ^ ^ ^ ^jg^ (^2„ _ ^2„) ^ ^jg^ (^1 _ ^2„+l)^ 

Sign (<^^-/ - xll-i') = sign (02'li - xl-_i) = sign (^^ _i " ) 

Sign(<^2n_,^2n) ^ sign (0^ - X^^+l ) 
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Now, since x""*"* € and since is satisfied on by definition, one obtains 

sign (0^+^+1 -a-^'+'+i) = sign ie{l---n}. (29) 

Combining tliis with the previous hst of cquahtics gives 

sign (0,1 - .t2"+i ) = sign (<^i ) , i e { 1 • • • n} , 

which exactly means a;^"+i = Tx^ G Z^, as expected. □ 

One may now consider the restriction of T to Z^. ActuaUy, the previous proposition imphes 
that chosing any intial condition £ W^, the trajectory {T™x | m G N} is contained in Z^, 
except maybe for x^. So, there is no loss of generality in chosing initial conditions in this 
invariant subset only. 

We can also verify that, on each wall all trajectories of system ([T7| are eventually contained in 
the region Z'. This fact is a consequence of proposition [5] and the following: 

Proposition 5. For all i (E {1 ■ ■ ■ 2n}, T°' {Z'-^) C Z\ 
Proof. Let x G Z\ which by definition is equivalent to 

Vj G { 1 • • • n} \ {.s,_ 1 } , sign (0;. - X, ) ^ sign (0} ) . (30) 

First, observe that in dimension ti = 2, = by definition and thus Zi(CT') = W^. There is 
nothing to prove in this case, and one thus assumes n 3 in the rest of the proof. 
From the alignment condition pT|) one gets, 

Yj^s., 0j-O7'x=(0}-i-x,)a,(x). 
Since aj{x) > 0, using ([501) leads to 

V? G {1 • • • n} \ 5j, sign (0} - 7fx) = sign (0}'i) = sign (0}). 

From (|24p it is immediate that 

sign (0;,_^ - G^^^x) = sign (</>;;:i^ )= sign (0:^_ J, 
which terminates the proof. □ 
This proposition allows us to consider only the restrictions: 

T'^' : Z*"i ^ Z' 

in the following. Figure [3] shows the domains Z* on a 3 dimensional example. 



4.3 Main result 

The main result of this paper will essentially consist in applying theorem [T] to feedback loop 
systems of the form (fTl]). or equivalently (fT7|) . Among the hypotheses to be checked, the 
concavity condition (C) will be verified with the aid of the following lemma. 

Lemma 1. Let p,q,r G R". Let also T : [0,p] —> [0, g] and M : [0,q] [0,r] be twice 
differ entiahle mappings satisfying condition (M) of theorem]^ Suppose that any of the second 
order derivatives of T and M are nonpositive, and that at least one of them is negative, for both 
mappings. 

Then, their composite MT : [0,p] ~-^ [0,r] satisfies (M) and (C). 
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3A 




Figure 3: The same walls as in figure [H Now, the invariant region [J - has been 
shaded in dark green. The reader may check on the figure that trajectories originating 
on any wall eventually fall into this region. 



Proof. The result is obtained from the ehain rule. Let x G [0,p]. Then Tx G [0, q] and 

D [MT) [x) = DM {Tx) DT{x). 

From this expression it is clear that the composite MT satisfies (M) if both T and M do. Now, 
given given any x < y in [0,p], condition {M) on M implies Tx < Ty, by monotonicity of each 
coordinate function T^, M;. 

Then the nonpositivity of second order derivatives gives DM{Ty) < DM{Tx). Then using both 
conditions (M) and (C) on T, and condition {M) on M , the chain rule expression leads to 

D {MT) {y) = DM {Ty) DT{y) < DM {Ty) DT{x) < DM {Tx) DT{x) = D {MT) {x). 

□ 

We are now in good position to state our main result. 
Theorem 2. Consider the map T : ^ defined in the previous section. 

• If n = 2, then \/x G Z^, T™x — *■ when m oo. 

• Ifn>2, then there exists a unique nonzero fixed point q = Tq. Moreover, q € Z^ and for 
every x G Z^ \ {0}, T™a; — > q as m ^ oo. 



Proof. The proof consists in checking the hypotheses of theorem [U which then yields precisely 
the above proposition. 

First, the upper corner of Z^ is the point (02^ . . . 0+), which belongs to R"~^. Second, it is clear 

from the expression ((23|) that each T° is C^. Hence, T is one the whole rectangle Zi. It 
follows in particular that DT{Q) is well defined. 

Since the proof is a little lengthy, it will be structured in three parts. 
Preliminary change of variables. 

Before checking the hypotheses of theorem [H it will prove convenient to apply a local trans- 
formation : Z^ ^ at each wall, such that J" C R"^^. This will allow us to deal with 
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nonnegative variables only. 

According to ([27|) . cr* denotes the sign of 0^, for j G Ji. From (j30p . 0^ and cc^ have the same 
sign, for G Z*. We also note cr] the sign of (p^j when j e {si-i, s.;}. Then, simply define the 
jth coordinate of this tranformation as 



(31) 



In other words, is just a multiplication by a diagonal matrix, whose jth entry is G {±1}- 
This implies that is invertible, and equals its inverse. Hence we also denote this inverse: 
: CP* — » Z*. One may also formulate J" explicitly: 



n 1^^- 



(32) 



It is tempting now to define maps M*-*-*, by the condition that the following diagram com- 
mutes, for any z G {1 • • • 2n}: 



(33) 













I' 







With the previous remark on this may also be written as M'-*' = p' T" p' ^. Composing the 



diagrams above for successive «, and since p^ 



id, one gets = Z^, and 



(34) 



For any ^ — [0,p] we write a;* = M'-*'x 
Each M(*) can be expressed coordinate-wise: 



^ 2. From proposition O one gets a;* G 



- (J'L ^ Xs 



3 7^ S i 



From their definition, the ajS clearly satisfy the alignment condition (pijl : cr* = ct* ^ for ^ Si 



It also comes from 



that: (tI 



Hence, one also has: 



where 



Monotonicity and concavity. 

One computes now: 



a," ^ X. 







.) j 7^ Si, 

10'.. I 



if i = fc 

if k ~ Si 
otherwise. 



(35) 



(36) 



and 



dxmdxk 



1] "i (f^l. ^ 



(^) = 



7s. 



7i 



10^ 



7s. ; (|0',J+a;,J 



2 



if j = fc, m = Si 
or fc = Si , m ~ j 

J if TO = A: = Si 
otherwise. 



(37) 
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One sees that M(*) is very similar to T*^ , so that equations above are strongly reminiscent of (fTU|) 
and pT|) . Let us evaluate the sign of these quantities. From x G and ([32]) it follows that 
10*1 —Xj ^ 0. It is clear also that 10*. | +Xsi > 0. It follows that, for all maps M^*-*, the nonzero 
terms of the jacobian are positive, while the nonzero derivatives of order 2 are negative, at any 
point X G T'^. 

Applied to the chain rule gives 



■DM'^^^x^] 



(38) 



Moreover, it seems from (j36p that the jacobians in the product above are the sum of a diagonal 
matrix and a matrix whose column number Si contains the only nonzero elements. 
Some care must be taken though, due to the particular indexing of matrices. Actually, from (j36p 
the diagonal entries of DM^'' are at first sight DJA^jj = aj(^al~^xl~^). However, the row 
subscripts of DM^'^ belong to {1 • • •ri}\{sj}, while the columns are indexed by {1 • • •n}\{si_i}. 

This implies that neither D3V[s.\si nor D3V[s._-^^si^i are defined. The exact shape of nonzero 
entries in Z?M^'-' depends on the order in which .Si_i and Si are sorted. In our case, table 
shows that either Si = Si-i + 1, or = 1 and Si_i = n, the latter occurring for i S {1, n + 1}. 
Let us represent the subscripts of nonzero entries in both cases. First, it i ^ {l,n+ 1}, both 
DM^*) and i:)M("+') have the following shape: 



(1,1) 



I 

(i,«0 



1) 

(Si_l, Si) 



: + l) 



{n,Si) {n,n) 

Now, ifiG{l,n+l} which implies that row 1 and column n are removed: 

(2,1) (2,2) 



(39) 



(n,l) 



(n - l,n - 1) 




(40) 



Since Si-i ranges over {1 • • • ti — 1} when i or n + i varies in {2 • • • ti}, it is not hard to deduce 
from ([5^ that both products: 

= DM(2")(a;2"-i)...i:>M("+2)(x"+i) and = M'-''\x''-^) ■ ■ ■ DM'-^\x^) 

have no zero entries, and are thus positive. Then, from (j40p one deduces that 

DT{x^) = DM^^'>{x^")U^ i:)M("+i)(a;")ni 

is positive as well. Hence we have just shown (M): 

yx^ey\ DT{x^)>Q. 

As for concavity, it follows from the negative terms in ((571) ^-^^ lemma [T] that 
D (m(*+i'm(')) {x'-^) < D (m(*+i)mW) (y'-i). 
Then a simple induction shows that DT. as expressed in satisfies condition (C): 
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Behaviour at the upper corner and the origin. 

Before studying the origin, let us show that Tp < p, as required by theorem [T] Here, p = 
{(f)2 ■ ■ ■ (f)n) is the upper corner of Z^, i.e. Z} = [0,p]. Now, let us consider the last maps 
composing T, as in the proof of proposition [H Namely, one deduces from that: 

Vx- e \ X- - ( 1 - 1 (x) ) , 

As usual, one identifies a-^"+^ and . Then for i e {2 ■ • • n + 1}, ai-i{x) £ (0, 1] gives: 

TtT2: = (l-a^-i(a;))0+i <0ti- 

Now, from equation ([2^ in the proof of proposition [H one knows that sign ((^"+* — Xsi_i) is 
conserved through the last steps of the cycle, so that 

Vi e {2 • • • n}, Vx G (t"^'?"'" ■ • • T'^"^') .t < 0+ 

Since in particular 

it follows that Tp < p. 



To terminate the proof, it remains to check whether the spectral radius of £'T(0) is greater 
than 1 or not. First, remark that T(0) = implies that each DM^*^ has to be evaluated at in 
expression (p8|) . Then from the expression of aj in (j36p . for z G {2 • • - n} the diagonal terms of 
each _DM'*-'(0) equal 1, except at — 1, i — 1). Here, it is given by the second line in just 
as the rest of column i — 1. The latter is given by: 



Vi e {2 • • • n} 



with the notation = Kj{a^). 
If the case n — 2. matrices I?M'-*'(0) arc in fact scalars of the form above. From equations ((T 
and ITTI) one obtains 



.2 _ 

— '^2 ^ — ' ' 1 — ' ^ 1 — ' 1 1 ^2 ~ — 9 — ^2 ' ctiiLi /\ — - 1% 

See figure m for the state space of two dimensional feedback loops. Then, gives 



and 



DT{0) 



Hir\ t\i-t Kin Hii-l 



tv-l tv<r) tx-t tvry 



1. 



Suppose now that n ^ 3. Before going further, let us remark the following fact: if a matrix A is 
such that A — Id > then the spectral radius p{A) > 1. This is a direct consequence of Perron- 
Frobenius theorem. Actually, p{A) is the maximal eigenvalue of A, or the maximal root of its 
characteristic polynomial, which we denote 7r^(Ar). One clearly has iTA-idi^X) = tta^X — 1), so 
that the roots of tta-IcI are exactly those of tta minus one. Applying Perron-Frobenius theorem 
to the positive matrix A — Id implies p{A) > 1. 

Hence to terminate the proof, it suffices to show that £'T(0) — Id > 0, or equivalently that 
the diagonal terms of DT(0) arc strictly greater than 1, since we have shown (M) already. We 
use from which it is possible to derive the general term of DT. 

The subscripts in and (jlUll arc those to be used in equation To compute on the 

other hand, it is more relevant to use the actual row and column numbers, as represented in (|39p 

(i) 

over the nonzero column, which is numbered Si-i. Hence in the expressions below, D3Vl -j^ refers 
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Figure 4: The state space of a feedback loop with only two variables, schematically. 
The notations are those explained in the text. 



to the entry j,k of matrix DM'-'mO), which is not necessarily (0). A simple induction 

dxk 

shows that 

n-l 

fel,...fe2re-l=l 

In particular, the diagonal terms can be conveniently denoted using the following abbreviation. 
Define the integer vector k = (fci . . . k2n-i) € {1 • ■ • n — and denote 



Note that the general form of terms in products above is the following: 

DmI^'^ . je{3---2n}. 

K2re+1- j fc2n + 2-j .J I ] 

Then, (^DT{x^)) equals the sum of products pj,. when k varies in the whole set {1 • ■ • n — 
Now, since each matrix DM*^*-* (0) is composed of zeros and positive terms, it follows that each 
nonzero product p]^ is positive. Hence, if for all i there is one k such that = 1, and at least 
one other nonzero pj^., the proof will be finished. 



To help intuition let us detail the shape of DM^'\0) and DM'^''+'\0) for n > i > 1, thanks 
to (011) and 

i - 1 
I 



For i e {1, n + 1}: 



r 1 



I 1 



(41) 



KM I 



(42) 
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Now, we claim that = 1, for 1 < i ^ 71 — 1, and 



(43) 



where k^ 



for 



1 — i ^ j ^ n, and 2n + l —i ^ j ^ 2n — 1, and kj = 1 otherwise. Actually, 



one has on the one hand: 



and on the other hand 



and 



DM 



(n+i+l) 



n+i+l 



^n+i+1 



and 



DM 



^2?i — i 5^' 



i+i 



(i) 



Moreover, for any j G {2?! + 1 — «---2n — 1}U{7t. + 1— it appears that -DM^ / = 1 

since i ^ Sj — I, and similarly, for j E {I ■ 
consequence, for k as in ([33]), 



- i} U {71 + 1 • • • 2n - i} one has DM[^\ = 1. As a 



i+l 1^1 



^i+1 



i+l 



K 



l+l 



Now we refer to the table (fT5|) of successive boxes, and recall that 
Q,i+i ^pj^QQ jg rather straightforward to check that 



(44) 

belongs to the domain 



sign( 



-sign ( 



for any i g {1 • • -n}. It follows immediately that sign(K"+') = — sign(K') as well. From this 
and the fact that any production term is fully determined by its sign, it follows that 



^i+l ^n+i+1 



„n+l 



= 1 



and 



^i+l ^-i+1 

.i+l ^.n+i+1 
"-i+l '^i+1 



1, 



whence Pk = 1 • 

Recall that the previous holds only for i > 1. Now, let k = (1 • • • 1). Then, 



and DM[1+^}1^ = ^^17'' 



^n+1 



while 



DM 



(n+2) 



^n+2 



,n+2 



Moreover for any other j, DM^-fl = 1, and thus 

A ^n+2 n 



and 



DM^^ 1 



dm'-II 



1^1 ^n+2 n+1 2 

A>Q rv-\ r\ir\ r\i-\ 



n+2 n+1 



1, 



Thus, we have found, for each i S {1 



(l---l,i- 



where kj = i for 2n + 1 



for the same reason as in ((H]). 

1}, a k such that = 1- Now, for any i > 1, let 
- i ^ j ^ 2n ~ 1. This leads to 

■ ■ DM['^^^ DM^J ■ ■ ■ DM^^}, 

and from ((^T|) and (|42p one checks readily that the above is positive: the entry (1, 1) is actually 
never zero in the first terms, the ith column of DM'*"*'^-' is nonzero, and the diagonal terms of 
the last terms in the product are positive as well. 

For i = 1, a positive pi. is found for instance with k = (2 • • ■ 2). Actually. (2, 2) entries arc never 
zero in (pij). and moreover DM\^1 and iJMg^J are also positive. □ 
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5 Conclusion 



The main result of this paper is theorem [21 which is a direct application of a theorem of 
Smith [13j , here refered to as theorem [TJ It is in fact very similar to the main theorem of 
Snoussi [11], since the latter states existence of a limit cycle where theorem [5] states existence 
and unicity, for the class of negative feedback loop systems of the form ([T7|) . The improvement is 
entirely obtained from the use of theorem[TJ This use is made possible by the natural occurrence 
of monotone, concave transformations in the context of piecewise afline gene network models of 
the form ([!]). In general, these transformations only appear on particular regions in phase space, 
as shown in section[3l Also, they are monotone and concave only up to a symmetry (or say with 
respect to the order induced by orthants which are not necessarily positive), which depends on 
the region on which they are defined. For negative feedback loop systems however, it is shown 
in section [mi that after a finite number of steps, all trajectories belong to one region, where the 
first-return dynamics is monotone and concave. It seems promising then, to study systems of a 
more general from than a negative feedback loop, using the same tools as in the present paper. 
For such systems, one should not concentrate on a single monotone concave transformation, 
but on several of them, and take into account the possible transitions between their domains of 
definition. 
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